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ABSTRACT 

Prestack  depth  migration  provides  a  powerful  tool  for  doing  velocity  analysis  in 
complex  media.  Both  prominent  approaches  to  velocity  analysis — depth-focusing 
analysis  and  residual-curvature  analysis,  rely  on  approximate  formulas  in  order 
to  estimate  velocity  errors.  Generally,  these  formulas  are  derived  under  the  as¬ 
sumptions  of  horizontal  reflector,  lateral  velocity  homogeneity,  or  small  oflTset. 
Therefore,  the  conventional  methods  for  updating  velocity  lack  sufficient  compu¬ 
tational  efficiency  when  velocity  has  large,  lateral  variations. 

Here,  using  a  perturbation  method,  we  derive  an  analytic  relationship  be¬ 
tween  residual  moveout  and  residual  velocity.  In  our  derivation,  we  impose  no 
limitation  on  offset,  dip,  or  velocity  distribution.  Based  on  this  formula,  we  revise 
the  residual-curvature-analysis  method  for  velocity  estimation  that  involves  recur¬ 
sion  and  iteration.  Furthermore,  this  formula  provides  the  s'^nsitivity  and  error 
estimation  for  migration-based  velocity  analysis,  which  is  helpful  in  explaining 
the  reliability  of  the  estimated  velocity. 


INTRODUCTION 

Prestack  depth  migration  that  can  handle  dipping  reflectors  and  lateral  velocity 
variations  is  robust  in  imaging  complex  structures.  In  order  to  process  data  by  this 
method,  one  often  needs  to  have  a  more  accurate  velocity  model  than  may  be  obtained 
from  simple  velocity-analysis  methods,  such  as  normal  moveout.  Meanwhile,  prestack 
depth  migration  itself  is  an  attractive  tool  for  doing  velocity  analysis  because  of  its 
high  sensitivity  to  the  velocity  model.  Two  approaches  to  migration  velocity  analysis 
have  been  developed:  depth-focusing  analysis  (DFA)  and  residual-curvature  analysis 
(RCA).  Depth-focusing  analysis  is  based  on  using  stacking  power  to  measure  velocity 
error.  Residual-curvature  analysis  is  based  on  residual  moveout  to  measure  velocity 
error. 

Migration  velocity  analysis  must  address  two  issues  to  succeed:  (1)  how  to  estab¬ 
lish  a  criterion  for  knowing  if  a  migration  velocity  is  acceptable;  (2)  how  to  update 
the  velocity,  if  it  is  unacceptable.  In  DFA,  a  migration  velocity  is  acceptable  if  the 
difference  between  migration  depth  and  focusing  depth  is  zero.  In  RCA.  a  migration 
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velocity  is  acceptable  if  the  difference  between  imaged  depths  from  different  offsets  is 
zero.  Quantitatively,  these  differences  can  be  used  to  update  the  velocity.  To  date, 
a  variety  of  updating  formulas  have  been  developed.  These  formulas  degrade  with 
increasing  complexity  of  media  (lateral  velocity  variation  and  reflector  dip),  because 
of  rough  approximations  used  to  estimate  velocity.  Although  iteration  generally  is 
helpful  in  obtaining  a  more  accurate  velocity,  a  good  approximation  not  only  reduces 
iteration  steps  but  guarantees  the  convergence.  The  conventional  formulas  for  updat¬ 
ing  velocity  were  derived  under  one  or  more  assumptions  as  follows; 

(i)  lateral  velocity  homogeneity; 

(ii)  small  oflFset; 

(iii)  horizontal  reflector. 

In  DFA,  all  three  assumptions  are  used  (MacKay  and  Abma,  in  1992).  In  RCA, 
the  assumptions  vary  with  migration  types.  For  common-shot  migration  or  common- 
receiver  migration,  .A.l-Yahya’s  formula  (1989)  used  all  three  assumptions.  Lee  and 
Zhang  (1992)  developed  a  formula  that  replaces  assumption  (iii)  by  the  assumption 
of  small  dip.  For  Common-offset  migration,  Deregowski’s  formula  (1990)  used  the 
first  two  conditions.  [This  was  proven  by  Liu  and  Bleistein  (1992).]  Both  DFA 
and  RCA  use  assumptions  (i)  and  (ii).  Under  the  assumption  of  a  small  offset,  the 
residual  moveout  can  be  approximated  by  a  hyperbola  or  a  parabola.  However,  this 
approximation  is  poor  when  the  velocity  has  an  obvious  lateral  variation. 

Conventional  approaches,  such  as  DFA  and  RCA,  inspect  information  on  each 
midpoint  individually.  This  may  not  be  reliable  when  conflicting  events  exist  because 
velocity  estimation  will  depend  on  which  event  is  used.  Furthermore,  the  conventional 
approaches  estimate  the  stacking  velocities  that  are  assumed  to  equal  the  rms  veloc¬ 
ities,  and  then  convert  the  rms  velocities  into  the  interval  velocities  by  Dix’  equation 
or  other  methods.  The  stacking  velocity  is  very  sensitive  to  the  lateral  variation  in 
velocity.  When  the  velocity  is  laterally  variant,  the  stacking  velocity  may  be  very 
different  from  the  rms  velocity. 

Because  of  limitations,  conventional  approaches  to  velocity  analysis  cannot  handle 
complex  media  (such  as  the  model  for  Marmousi  data).  To  overcome  these  limitations, 
some  geophysicists  used  a  macro-model  method  for  estimating  the  interval  velocity 
directly  .  A  macro-model  consists  of  velocities  and  velocity  interfaces.  An  initial 
macro-model  is  determined  from  conventional  techniques;  the  max:ro-model  is  updated 
iteratively  by  using  depth-focusing,  residual-curvature  analysis,  or  both.  The  whole 
process  may  be  accomplished  with  interpretation  based  on  geological  knowledge.  The 
present  macro-model  method  has  two  significant  drawbacks.  The  first  one  is  that 
velocity  distribution  is  assumed  to  be  constant  in  layers.  In  fact,  velocity  may  vary 
within  a  block.  The  second  one  is  that  no  efficient  formula  is  used  to  update  velocity. 
The  criterion  used  now  for  updating  velocity  is  only  qualitative,  which  results  in 
inaccuracy  and  inefficiency  in  velocity  estimates. 

In  this  paper,  we  present  a  new  method  for  velocity  analysis  that  revises  the  RCA 
method.  Based  on  perturbation  theory,  we  derive  a  quantitative  relationship  between 
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residual  moveout  and  error  of  velocity  that  is  valid  for  any  offset,  dip  and  velocity 
distribution.  A  similar  formula  was  proposed  by  Lafond  et  al.  (1993).  However, 
that  formula  is  restricted  to  const  ant- velocity  layers.  In  contrast,  our  formula  can 
handle  any  parameterization  of  velocity.  To  apply  this  formula,  it  is  necessary  to 
calculate  some  quantities  that  depend  on  the  specular  source-receiver  pairs  for  each 
image  point.  Here,  we  use  multiple  outputs  of  Kirchhoff  migration  to  determine  these 
(juantities  efficiently. 

.\lthough  prestack  migration  has  superior  advantages  to  deal  with  complex  me¬ 
dia,  its  weaknesses  should  be  considered  in  application.  Besides  the  costly  computa¬ 
tion,  complexity  of  medium  may  degrade  the  effectiveness  of  the  prestack  migration 
method.  Here,  we  discuss  the  sensitivity  of  migration-based  velocity  analysis,  and 
show  what  factors  affect  the  sensitivity,  so  that  we  may  estimate  the  velocity  error 
involved  in  the  velocity  analysis.  Furthermore,  this  sensitivity  analysis  provides  an 
indicator  of  when  a  priori  geological  knowledge  is  necessary  for  estimating  reasonable 
velocities. 


VELOCITY  ESTIMATION 


If  an  incorrect  velocity  is  used  in  prestack  migration,  the  imaged  depths  from  dif¬ 
ferent  offsets  at  a  common-image-gather  will  differ  from  each  other.  In  this  situation, 
a  residual  moveout  is  observed  in  migrated  data.  The  principle  of  velocity  analysis 
is  to  correct  the  velocity  so  that  the  imaged  depths  at  each  common-image-gather 
are  close  to  for  each  other.  For  this  purpose,  one  iteeds  a  quantitative  relationship 
between  the  residual  moveout  and  the  velocity  error.  Here,  we  use  the  perturbation 
method  to  derive  the  representation  for  the  residual  moveout. 

We  denote  by  A'  a  2-D  vector,  A'  =  (x,  z)  .  Let  x,  be  the  source  position  and  Xr 
be  the  receiver  position  on  the  surface.  For  any  point  A'  below  the  surface,  r3(x3,  A') 
or  Tr(X,Xr),  respectively,  denote  traveltimes  from  x,  to  A',  or  A”  to  x^. 

Suppose  we  know  the  total  reflection  travetime  function  T{y,  h)  (therefore,  dTjdy) 
that  depends  on  midpoint  y  and  half-offset  h.  Given  a  velocity  function  v{x,  c),  then, 
for  each  h.  the  reflector  is  determined  by 


r3(x3,  A')  -I-  r^A'.x,.)  =  T(y,h), 


dr.  dVr 


-I- 


(1) 

(2) 


dy  '  dy  dy' 

where  A'  =  (x,  c)  is  the  point  on  the  reflector.  The  reflector  is  related  to  offset;  for  a 
fixed  image  location  x,  the  imaged  depth  z  can  be  computed  that  is  a  function  of  h. 
If  i;(x,  z)  equals  the  true  velocity,  then  the  imaged  depth  2  is  independent  of  offset  h; 
otherwise — for  incorrect  velocity — ^  varies  with  offset  h.  Consequently,  the  imaged 
depth  provides  information  on  velocity. 

Suppose  that  the  velocity  distribution  v  is  characterized  by  a  parameter  A. 

V  =  v(X;X). 
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For  example,  when  v(x,  z)  =  vq  a  ax  +  bz,  A  is  either  I'o,  a,  or  b. 

For  a  fixed  image  location  x,  we  differentiate  equation  (1)  with  respect  to  A. 
Noticing  that  y  and  are  functions  of  A,  then 


dr,  ^ 

dy  dy  dX 


dr,  dxr 
'dX^^ 


+ 


dr,  dxr  dz 
dz  dz  dX 


dT  dy 
dy  dX 


(3) 


By  using  equation  (2),  the  first  term  of  the  left  side  in  equation  (2)  is  balanced  by 
the  right-hand  term.  Therefore, 


dr,  dxr  dz 
dX 


dr,  dTr 
dX  dX' 


From 


dr, 

dz 


cosO, 


dtr 

dz 


cos  9r 


V  az  V 

where  9,  or  9^  are  angels  between  the  raypath  from  the  source  or  the  receiver,  and 
the  vertical  at  A',  we  have  that 


That  is, 

where, 


cos  9,  +  cos  9,  dz  dr,  dr^ 

V  dX  dX  dX 


dz 

dX 


9{x,h), 


g{x,h)  =  - 


dr,  dTr 
dX  dX 


f(A;X) 
cos  9,  4-  cos  9r 


(4) 

(5) 


Suppose  that  the  true  parameter  is  A*  and  the  true  depth  is  z'.  If  there  is  a 
small  perturbation  AA  =  A*  —  A  between  the  true  parameter  and  the  parameter 
used  in  migration,  then  the  imaged  depth  will  have  a  corresponding  perturbation 
Az{x,h)  fn  z*  —  z{x,h).  From  equation  (4),  Az  can  be  represented  by 


Az  =  g{x,h)AX. 


(6) 


Equation  (6)  is  true  for  an  arbitrary  velocity  distribution,  any  dipping  reflector, 
and  any  offset,  which  is  significantly  different  from  the  limited  result  of  conventional 
RCA. 

When  the  velocity  distribution  is  characterized  by  multiple  parameters.  A  =  (Aj, 
A2,  ...,  A„)^,  the  imaged-depth  perturbation  will  be  affected  by  the  perturbations  of 
all  these  parameters.  Therefore,  equation  (6)  is  modified  by 

A.-(T,ft)  =  Y,  =  £5,(i,*)AA„  (7) 

i=l  ,=  1 
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where 

r(A:X) 
cos  +  cos  dr 

If  we  could  solve  for  AA.  the  true  parameters  can  be  estimated  by 

A*  =  A  +  AA. 

Notice  that  the  left  side  of  equation  (7)  involves  the  unknown  so  it  is  necessary 
to  remove  in  this  equation.  The  true  depth  can  be  approximately  replaced  by  the 
corrected  imaged  depth  +  A~, 


f/.i-r  Ji)  =  - 


pTr 

dX,  OX, 


c’  «  ;(x,  h)  +  Ar(j‘,  h)  =  z{x,  h)  +  ^  g,{x,h)AXi. 

1=1 

Because  z’  is  independent  of  offset,  the  corrected  imaged-depths  from  different  offsets 
should  be  close  to  for  each  other.  Suppose  that  there  are  offsets  /q,  /12,  ...,  and 
image  locations  jj,  •••<  -Tas  f^ien 


+AiW=.-f' +5:»!fAA,, 


1=1 


where 

=  z{Xk,hj), 

=  Az{xk,hj), 

9ii^  =  gi{xk,hj). 

We  seek  AA,’s  such  that  the  corrected  imaged  depths  have  the  minimum  variance; 

K  m  _ 

^  H  =  min,  (9) 

fc=i j=i 

iW  _  J>^)  Ak))T 

~  —  V^l  5  ~2  5  "-I  "m  )  5 

Ai‘*^^  =  (Acl*^\Acf,...,AcW)^. 

Here  we  use  the  overline  to  denote  the  mean  value  of  a  vector  over  the  offset  index. 
For  example, 


i.e.. 


where 


1  "• 


{k) 


Indroduce  the  matrix  and  vector, 


h- 


o 
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where 


then,  as  shown  in  Appendix  A.  the  solution  of  equation  (9)  must  satisfy  the  linear 
equation, 


K 


K 


(10) 


LJb=i 


;fc=l 


Specifically,  if  there  is  only  one  parameter  to  be  determined,  equation  (10)  will  have 
an  explicit  solution 


where 

and 


=g{xk,hj), 


(11) 


If  the  corrected  imaged-depths  are  not  close  enough  to  each  other,  we  implement 
iteration  to  obtain  more  accurate  parameters.  The  iteration  stops  when  the  variance 
achieves  a  given  accuracy. 


Calculation  of  the  function  gi 

The  function  gi{y,h)  involves  the  derivatives  of  traveltimes  with  respect  to  the 
parameter  A.  In  the  eikonal  equation. 


1 

v^{x,  z) ’ 


if  we  take  the  derivative  with  respect  A,  then 


where 


dyidr  dydr  _  \  d  (  1  \ 

dx  dx  ^  dz  dz  2  0A  z)  J  ' 

dr 


(12) 


For  each  source  or  receiver,  drjdX  is  determined  by  solving  equation  (12).  Therefore, 
given  an  imaged  point  {x,z)  and  a  specular  source-receiver  pair  Xg  and  Xr,  we  can 
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calculate  g,  from  formula  (8).  However,  there  is  not  an  explicit  formula  to  represent 
the  specular  source-receiver  pair  from  the  imaged  point  for  a  complex  medium.  To 
solve  this  problem,  here,  we  use  the  Kirchhoff  integral  to  calculate  g,.  In  the  Kirchhoff 
summation,  we  calculate  two  inversion  outputs  which  have  two  different  amplitudes. 
One  is  the  original  amplitude:  the  other  is  the  original  one  multiplied  by  the  quantity 

Thus,  the  ratio  of  the  amplitudes  of  the  two  outputs  will  give  at  the  specular 
source-receiver  position  according  to  the  stationary-phase  principle,  without  requiring 
knowledge  of  the  specular  source-receiver  pair.  This  is  the  same  technique  as  was  used 
to  determine  the  angle  of  reflection  in  Kirchhoff  inversion.  [See  Bleisten  et  al.  (1987).] 

Parameterization  of  Velocity 

Although  equation  (7)  holds  for  any  velocity  distribution,  the  solution  will  be 
underdetermined  if  too  many  unknown  parameters  are  involved.  Consequently,  it  is 
essential  to  characterize  the  velocity  distribution  by  choosing  appropriate  parameters. 
Conventionally,  one  assumes  that  a  velocity  model  consists  of  the  construction  of  the 
macro-model  (constant  velocities  and  velocity  interfaces).  The  interfaces  divide  the 
whole  model  into  a  number  of  blocks.  Here,  we  replace  constant  velocity  in  one  block 
by  a  linear  function  that  is  characterized  by  three  parameters: 


Al  -I-  A2(2  —  ^o)  +  '^3(^  ~  -To)) 

where  (iq,  ^o)  is  a  reference  point.  Thus,  the  velocity  distribution  is  written  in  a  form 
of 

v(x, :)  =  z'o(x, .:)  -f  Al  +  A2(2  -  2ro)  +  A3(x  -  Xo),  (13) 

where  vq  is  a  background  velocity. 

An  imaged  depth  only  depends  on  the  velocity  above  it,  except  for  turning  rays. 
Therefore,  a  recursive  algorithm  is  possible  to  determine  velocity  in  an  individual 
block.  We  start  from  the  block  nearest  surface.  In  each  block,  iteration  is  used  to 
calculate  velocity  parameters.  Given  an  initial  guess  for  A^’s,  common-offset  depth 
migration  is  implemented  to  obtain  imaged  depths  and  g,(x,/i)  in  equation  (8)  for 
each  midpoint  and  each  offset.  Equation  (10)  will  give  a  correction  of  the  parameters. 
Then  by  using  the  updated  parameters  as  initial  guess,  we  correct  the  velocity  again 
until  convergence  is  achieved.  After  velocity  analysis  in  one  block,  we  migrate  data 
with  the  correct  velocity,  and  then  pick  the  block  interface  from  the  imaged  structure. 
When  we  finish  determining  velocity  and  velocity  interface  in  one  block,  we  will  repeat 
the  same  procedure  to  the  next  block. 

SENSITIVITY  OF  MIGRATION  VELOCITY  ANALYSIS 

Velocity  analysis  by  prestack  migration  uses  the  difference  between  the  imaged 
depths  from  different  offsets  to  correct  the  velocities,  which  is  represented  by  equa¬ 
tion  (10).  If  the  variance  defined  by  (9)  is  zero,  we  conclude  that  the  velocity  is 
correct.  However,  one  cannot  obtain  exactly  zero  variance.  Many  factors  result  in 
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a  non-zero  variance;  noise  in  the  input  data,  nonacoustic  properties,  and  inaccurate 
description  of  velocity  distribution,  so  on.  Even  the  variance  is  apparently  zero,  it 
is  actually  not  zero  because  of  the  errors  in  picking  the  imaged  depths.  The  imaged 
depths  in  the  variance  are  picked  on  the  migration  output,  so  that  the  position  error 
in  imaged  depths  is  controlled  by  the  resolution  of  the  migration  output,  which,  in 
turn,  depends  on  wavelength,  among  other  things. 

Quantitatively,  the  matrices  in  equation  (10)  that  depend  on  the  functions,  5,, 
can  be  used  to  describe  the  sensitivity  of  the  velocity  error  AA  to  the  variance  error. 
Here,  we  will  derive  analytical  representations  for  the  simplest  cases.  For  simplicity, 
we  assume  that  the  velocity  v{x,  r)  consists  of  a  constant  background  velocity  t’o  and 
a  perturbation  that  is  a  linear  function  of  depth  in  one  block.  Moreover,  w^e  assume 
that  the  upper  boundary  of  the  block  is  a  horizontal  line,  c  =  d;  i.c., 

v{x,  z)  =  vq  +  o{Ai  +  \.i{z  -  ^o)),  (14) 

where  a  =  0  for  z  <  d  and  a  =  1  for  z  >  d.  The  initial  guesses  are 

Aj  =  A2  =  0. 


Suppose  that  the  true  parameters  are  Aj  and  A2,  and  the  reflector  is  a  horizontal 
segment  at  a  depth  2*.  As  shown  in  Appendix  B,  we  obtain  the  representations  for 
the  Qi's  in  equation  (8) 


gi{x,h) 


A  z^  z  —  d 


(15) 


g2{x,h) 


/l^-F2^(z-2o)^-(d-2o)2 

2vo 


(16) 


From  the  above  two  equations,  we  conclude  that 


g2{x,h) 


{z  -  zq)^  -  {d-zof 
2{z  ~  d) 


gi{x,h). 


This  means  that  the  coefficients  in  equation  (8)  are  proportional  for  all  image  locations 
and  all  offsets,  so  that  we  cannot  solve  for  Aj  and  A2  separately.  Thus,  we  have  a 
conclusion:  The  parameters  Ai  and  X2  cannot  be  determined  at  the  same  time  if  only 
one  reflector  segment  is  used  in  the  velocity-analysis  process. 

This  conclusion  shows  that  one  should  avoid  solving  for  Ai  and  A2  at  the  same 
time  unless  well-separated  reflector  segments  are  used  for  velocity  analysis. 

If  the  parameter  A2  is  given  (i.e.,  AA2  =  0),  then  equation  (8)  is  simplified  to 

/  \  2  —  d 

Az{xji)  =  +  Ij  — — AAi.  (17) 


When  AAi  is  small,  z{x,h)  «  2’.  Therefore,  for  any  different  two  half-offsets  /i2  and 
hi  (/i2  >  hi),  the  difference  of  imaged  depths  from  these  two  offsets  are 


z{x,h2)  -  z{x,hi) 


h\ 


h? 


[(2*)-^  (2-)2J  CO 


d 


AAi. 
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i.e.. 


AA,. 


(18) 


Equation  (18)  shows  that  the  error  in  Aj  is  proportional  to  the  thickness  of  the  block 
layer,  which  is  consistent  with  results  in  the  paper  of  Liu  and  Bleistein  ( 1990).  When 
there  is  a  thin  layer,  the  velocity  estimation  will  become  unreliable,  so  that  migration- 
based  velocity  analysis  cannot  handle  thin  layers. 


Similarly,  if  the  parameter  \i  is  given  (i.e..  AA]  =0).  then 


6z{x,hi,h-2)  = 


(hl-hi)iz'  ~d){z'  Ad-2zo) 
2(r-)2eo 


AA.,. 


(19) 


Equation  (19)  shows  that  the  error  in  A2  is  determined  not  only  by  the  thickness  of 
the  block  layer  but  by  the  difference  between  the  central  depth  of  the  block  layer  and 
the  reference  depth  Cq. 

We  have  two  kinds  of  the  depth  errors;  Ac  and  6z.  The  former  is  the  difference 
between  the  true  depth  and  the  imaged  depth,  which  reflects  migration  error  due  to 
the  velocity  error:  the  latter  is  the  difference  between  the  imaged  depths  from  different 
offsets,  which  results  in  the  velocity  error  in  velocity  analysis.  The  ratio  of  these  two 
depth  errors  is  given  by  equations  (17)  and  (18), 


Az{x,h)  _/i2-|-(c’)^ 
^  ~  <5c(x./li,/l2)  (/i2“M) 


(20) 


The  ratio  7  indicates  the  migration  error  due  to  the  resolution  in  migration  output. 
Equation  (20)  shows  that  7  increases  as  the  depth  increases. 

For  a  general  case,  we  do  not  have  analytical  representations  for  the  g.'s.  But  we 
can  calculate  the  Qi's  numerically  from  which  to  estimate  the  velocity  error.  In  fact, 
if  there  is  only  one  parameter,  equation  (11)  gives  an  explicit  error  estimation. 


|AA|  < 


(21) 


where  the  Cauchy-Schwarz  inequality  is  used. 

COMPUTER  IMPLEMENTATION 

We  applied  our  velocity  analysis  technique  to  synthetic  data.  The  prestack  migra¬ 
tion  is  implemented  by  Liu’s  program  (1993).  The  velocity  model  shown  in  Figure  1 
consists  of  linear  velocity  functions  and  interfaces.  Synthetic  seismic  traces  are  gen¬ 
erated  from  this  model,  with  5  offsets  ranging  from  100  meters  to  900  meters.  Two 
of  the  common  offset  gathers  are  shown  in  Figure  3. 
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The  velocity  analysis  process  is  outlined  as  follows; 

In  the  first  layer,  the  initial  guess  of  velocity  is  1500  ni/s.  With  two  iterations, 
the  velocity  is  2017  in/s.  The  true  value  is  2000  m/s. 

In  the  second  layer,  we  assume  that  velocity  has  the  foi  u  v(z)  =  2017-t  A(c  — 500) 
m/s,  which  guarantees  a  continuous  velocity  across  the  first  interface.  The  initial 
guess  is  A  =  0.  With  two  iterations,  the  velocity  is 

c(c)  =  2017  +  1.035(c  -  500)  =  1500  +  1.035r  m/s. 

The  true  velocity  is  t;(c)  =  1500  +  l.Or  m/s. 

In  the  left  third  layer,  we  assume  that  velocity  has  the  form  r{T,c)  =  1000  + 
1.5c  +  Xx  m/s.  The  initial  guess  is  A  =  0.  With  one  iteration,  the  velocity  is 
v(x,z)  =  1000+  1.5c  —  O.Ij"  m/s.  The  true  velocity  is  e(T,c)  =  1000+  1.5c  —  0.2x 
m/s.  Despite  the  error  in  lateral  variation,  the  stopping  criterion  is  met  since  the 
difference  of  the  imaged  depths  is  apparently  zero.  Therefore,  iteration  ceases. 

In  the  right  third  layer,  we  assume  that  velocity  has  the  form  i'(c)  =  2535 +  A(c  — 
1000)  m/s.  The  initial  guess  is  A  =  1.C35  With  one  iteration,  the  velocity  is 

c(c)  =  2535  +  1.55(c  -  1000)  =  985  +  1.55c  m/s. 


The  true  velocity  is  v{z)  =  1000  +  1.5c  m/s. 

In  the  fourth  layer,  the  velocity  is  a  constant.  The  initial  guess  is  u  =  4000.  With 
one  iteration,  the  velocity  is  3586  m/s.  The  true  velocity  is  u  =  3500  m/s. 

•After  velocity  analysis,  we  obtain  the  velocity  model  shown  in  Figure  2.  With  this 
velocity  model,  we  implement  migration,  and  the  result  is  shown  in  Figure  5  which  is 
close  to  the  migration  result  with  the  true  velocity,  shown  in  Figure  4. 

.Also,  we  test  the  computation  of  gi  in  the  first  layer.  The  numerical  is  computed 
from  the  ratio  of  the  amplitudes  in  the  KirchhofF  migration  outputs.  Then  we  compare 
the  numerical  gi  with  the  true  value  calculated  in  equation  (17).  In  our  test,  I’o  =  500 
m/s  and  d  =  0.  The  result  listed  in  Table  1  shows  that  our  numerical  \'alue  of  gi  is 
accurate  for  all  offsets. 


Table  1.  Test  for  <7i.  The  unit  of  offset  and  depth  is  m 


Offset 

Imaged  depth 

Numerical  gi 

True  gi 

100 

370 

0.249 

0.251 

300 

360 

0.278 

0.282 

500 

330 

0.345 

0.346 

700 

290 

0.476 

0.475 

900 

210 

0.773 

0.783 

/ 


s. 
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CONCLUSION 

Velocity  analysis  by  prestack  migration  can  liandle  complex  media  for  which 
conventional  approaches  lack  the  efficiency  in  updating  velocity.  The  perturbation 
method  is  used  here  to  derive  an  analytical  relationship  between  the  residual  moveout 
and  the  error  in  the  background  velocity.  This  method  is  not  limited  by  offset,  dip, 
or  velocity  distribution.  Based  on  that,  we  can  estimate  velocity  elhcientlv.  Further¬ 
more.  this  relationship  helps  us  to  estimate  the  velocity  error  from  the  data  resolution. 
This  velocity  analysis  technique  also  can  bo  generalized  to  3D  and  to  converted  waves. 
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Fig.  3.  Synthetic  data;  (a)  with  offset  of  100  meters  and  (b)  with  offset  of  900 

meters. 
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Fig.  4.  Prestack  migration  with  the  true  velocity  model  in  Figure  1.  The  offset  is 

100  meters. 
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Fig.  5.  Prestack  migration  with  the  velocity  model  in  Figure  2.  The  offset  is  100 

meters. 
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DERIVATION  OF  EQUATIONS  (10) 

By  using  the  matrix  and  vector  notations  on  page  5  and  6.  we  have 


1 

1  m  n 


J=lt=l 


n  /  1  ^ 

=  E  AA,iE«!r 


"■“l 


=  t  (aA,9,!‘')  ,  (22) 

and 

+  A2j*^^  —  Ai*^* 

=  .f  +  -AA,  (s!-’  - iP)  ■  (23) 

i=l  ^  ^ 

Define  a  function 

K  m  ■  2 

/(AA)  =  E  E  (^J  +  -  i"-'  +  'Ail*')  . 

fc=l j=l 

Thus,  finding  the  minimum  of  /(AA)  is  equivalent  to  that  its  gradient  with  respect 
to  AA  equal  6,  i.e., 

/  =  l,2,...,n.  (24) 

aAA/ 

By  using  equation  (23), 

/(AA)  =  E  E  k’  -  ^ + E  AA'  (s!f  -  iP)|  ■ 

fc=ii=iL  1=1  ^ 


=  2  E  E  U'"”  -  +  E  AAi  (dP  -  iP)]  (sif  -  Si' 

it=ij=iL  1=1  ^  ^ 

=  2EE(.p-iJ^(sP-^ 


fc=l  j=l 


+2 E  E E  AA,  (jif  -  9l‘’)  (sii’  -  j)‘')  • 

Jk=lj=li=l  ^  ^  ^ 
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The  conditions  (24)  imply  that 


K  m  n 


E  E  E  (4’  -  (slf  -  iF)  =  -  E  E  (--F  -  (sF  -  9.“’) : 

k=ij=ii=\  ^  ^  ^  '  k=ij=i  ^  ' 

i.e., 

=  (26) 
.*:=1  k=l 

This  completes  the  verification  that  the  solution  of  equation  (10)  provides  the 
minimum  of  the  left  side  in  (9). 

DERIVATION  OF  EQUATIONS  (14)  AND  (15) 


The  solution  of  equation  (12)  can  be  represented  by 


,dv{x,z) 


where  L  is  the  raypath.  Here,  to  simplify  notations,  we  suppose  that  all  derivatives 
with  respect  to  A  are  evaluated  at  Ai  =  A2  =  0.  When  the  background  velocity  is  a 
constant  vq  as  in  (13),  equation  (27)  becomes 

dr  1  f‘  _2dv(a;,  2')  ,  , 

m  = - n  /  ^0  V  ^  28 

aX  cos  0  Jo  dX 

where  0  is  the  angel  between  the  raypath  and  the  vertical.  For  the  velocity  function 
in  equation  (14),  if  2  >  d,  then 

dv(x,z') 


if  2  <  d,  then 


dvjx,  z')  _  j 
0X2  “ " 

dv{x,z')  dv{x,z') 


Therefore, 


dAi  cosd  Jd  Vq 

dr  ^ _ ]_  n  z'  -  zq^j  _ _ 1  (2  -  20)^  -  (d  -  20)^ 

dXo  cos  0  Jd  t’n  COS0  2v^ 


1  z  —  d 
cos  0  Vq 


Id  eg 
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Because  the  reflector  is  horizontal,  the  raypaths  from  the  source  and  the  receiver  are 
symmetric.  We  have 

cos6.  =  cos6r  =  z. 


Substituting  the  above  formula  and  equations  (29)  and  (30)  into  equation  (8),  we 
obtain 

z  —  d 

giix.h)  = - - - , 


+  z- {z  -  zq)^  -  {d  -  zo)"^ 
9,ix,h)  =  — - — - . 

This  completes  the  verifications  of  equations  (15)  and  (16). 
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